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Abstract 

We make a global vacuum neutrino oscillation analysis of solar neutrino data, in- 
cluding the seasonal and energy dependence of the recent Super-Kamiokande 708- 
day results. The best fit parameters for v e oscillations to an active neutrino are 
5m? = 4.42 x 10~ 10 eV 2 , sin 2 29 = 0.93. The allowed mixing angle region is con- 
sistent with bi-maximal mixing of three neutrinos. Oscillations to a sterile neutrino 
are disfavored. Allowing an enhanced hep neutrino flux does not significantly alter the 
oscillation parameters. 



1 Introduction 



For quite some time, measurements of solar neutrinos jl], have indicated a suppression 
compared to the expectations of the standard solar model (SSM) J|. This suppression may 
be explained by assuming that v e from the Sun undergo vacuum oscillations as they travel 
to the Earth ||, |]. Recent data from the Super-Kamiokande (SuperK) experiment || also 
exhibit a seasonal variation above that expected from the l/(distance) 2 dependence of the 
neutrino flux, which if verified would be a clear signal of vacuum oscillations 0, |3|. 

In this Letter we determine the best fit vacuum oscillation parameters to the combined 
solar neutrino data, including the 708-day SuperK observations. For oscillations to an active 
neutrino there are subtle changes in the allowed regions compared to fits with earlier SuperK 



data J9], [L0| . Oscillations to a sterile neutrino are disfavored. We examine the possibility 
that the solar hep neutrino flux is enhanced compared to the SSM, and find that the best 
fits are only marginally changed. We also find a solution with very low 5m 2 ~ 6 x 10~ 12 eV 2 
for oscillations to a sterile neutrino or with an hep enhancement. 



2 Fitting procedure 

The fitting procedure is described in detail in Ref. [[|. The solar data used in the current 
fit are the 37 CI (1 data point) and 71 Ga (2 data points) capture rates Q, the latest SuperK 
electron recoil energy spectrum with E e in the range 5.5 to 20 MeV |J (18 data points), 
and the latest seasonal variation data from SuperK for E e in the range 11.5 to 20 MeV 
(8 data points). For v e oscillations to an active neutrino [y^ or v T ) we take into account 
the neutral current interactions in the SuperK experiment. We fold the SuperK electron 
energy resolution || in the oscillation predictions for the E e distribution. For all rates that 
are annual averages we integrate over the variation in the Earth-Sun distance. We do not 
include the SuperK day/night ratio, since it is unity for vacuum oscillations. 

To allow for uncertainty in the SSM prediction of the 8 B neutrino flux, we include as a 
free parameter f3, the 8 B neutrino flux normalization relative to SSM prediction. In our fits 
with non-standard hep neutrinos we also include an arbitrary hep normalization constant, 
7. For the SSM predictions we adopt the results of Ref. ||. 

3 Solutions with oscillations to an active neutrino 

Figure la shows the 95% C.L. allowed regions for the combined radiochemical ( 37 C1 and 
71 Ga) capture rates and the SuperK electron recoil spectrum, including its normalization. 
For the 95% C.L. region in a three-parameter fit (5m 2 , sin 2 29, (3) we include all solutions 
with x 2 < Xmin + 8- On the boundary curves of the fit to the radiochemical data, the 
faster oscillations are due to the lower energy neutrinos (mainly pp and 7 Be) , and the slower 
oscillations are due to the 8 B neutrinos. The five regions allowed by the SuperK spectrum 
data correspond roughly to having a mean Earth-Sun distance equal to (in increasing order 
of 5m 2 ) ~, §, |, ~, or | oscillation wavelengths for a typical 8 B neutrino energy. Following 
the notation of Ref. [0], we label these regions A, B, C, D, and E, in order of ascending 5m 2 . 
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In Figure lb we show the allowed regions for the combined radiochemical and SuperK 
spectrum data (the solid curve). Also shown in Fig. lb is the region excluded by the seasonal 
SuperK data at 68% C.L. (we show the 68% C.L. region since almost none of the parameter 
space is excluded by the seasonal variation at 95% C.L.). 

Finally, in Fig. lc we show the 95% C.L. allowed regions obtained with all of the data 
(radiochemical, SuperK spectrum, and SuperK seasonal). Only regions A, C, and D are 
allowed at 95% C.L. by the combined data set. The best fits in each of the subregions are 
shown in Table [I]. Region B, which was allowed in previous fits [^, 10|, is now excluded. The 
overall best fit parameters are in region C 

5m 2 = 4.42 x 1(T 10 eV 2 , sin 2 29 = 0.93, (5 = 0.78 , (1) 

with x 2 1 DOF = 33.8/26, which corresponds to a goodness-of-fit of 14% (the goodness-of- 
fit is the probability that a random repeat of the given experiment would observe a greater 
X 2 , assuming the model is correct). The next best fit is in region D with 

6m 2 = 6.44 x 10" 10 eV 2 , sin 2 29 = 1.00, = 0.80, (2) 

with x 2 /DOF = 36.7/26, which corresponds to a goodness-of-fit of 8%. Finally, region A 
is also allowed, with best fit 

5m 2 = 6.5 x 10" 11 eV 2 , sin 2 29 = 0.70, = 0.94, (3) 

with x 2 /DOF = 38.4/26, which corresponds to a goodness-of-fit of 6%. Regions C and D 
are consistent with bi-maximal flTTf or nearly bi-maximal [ E^] three-neutrino mixing models 



that can describe both solar and atmospheric neutrino data. Regions B and E are excluded 
at 95% C.L. 

We see from Table [I] that generally the higher 5m 2 solutions fit the E e spectrum better. 
This is evident in Fig. 2, which shows the E e spectrum predictions for the solutions in Table [TJ, 
and the latest SuperK data |J. The higher 5m 2 solutions do worse for the radiochemical 
experiments, primarily because the seasonal variation in the oscillation argument is larger in 
these cases for the 7 Be neutrinos, causing more smearing of the oscillation probability. The 
annual-averaged suppression of 7 Be neutrinos is about 93% for solution A, 66% for solutions 
B and C, and 57% for solutions D and E. The best fit in each region generally lies near a local 
minimum for the 7 Be neutrino contribution, as might be expected from model independent 
analyses without oscillations that indicate suppression of 7 Be contributions [|T^]. 

Results for the seasonal variation of the solutions in Table [I] are shown in Fig. 3, along 
with the latest SuperK data 11. Also shown is the best fit for no oscillations with arbitrary 
8 B neutrino normalization; the latter has a small seasonal dependence due to the variation 
of the Earth-Sun distance. Current seasonal data do not provide a strong constraint, but 
clearly solutions A, C, and D fit the data better than the no-oscillation curve. 

Also included in Table [l] is a solution with very low 5m 2 , which we call solution Z, that 
was pointed out in Ref. [|TJ]]. It is excluded at 95% C.L., but marginally allowed at 99% C.L. 
(as is solution B). It corresponds roughly to a mean Earth-Sun distance equal to | of the 
oscillation wavelength (maximal supression) for 7 Be neutrinos. Although it does very well 
in suppressing the 7 Be neutrinos (by about 94% after accounting for seasonal averaging), 
it does not show a significant seasonal effect beyond that provided by the variation of the 
Earth-Sun distance (see Fig. 3b). 
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4 Solutions with oscillations to a sterile neutrino 

A similar analysis may be made with solar v e oscillating into a sterile neutrino species. Since 
the sterile neutrino does not interact in any of the detectors, it is harder to reconcile the 
differing rates of the 37 C1 and SuperK experiments ]7j . The best fit parameters in each of the 
regions are shown in Table 0. For sterile neutrinos the overall best fit parameter is again in 
region C with x 2 /DOF = 42.4/26, which is excluded at 97.7% C.L. Therefore oscillations to 
sterile neutrinos are highly disfavored. Interestingly, the fit for solution Z is comparable to 
the other best-fit solutions in the sterile case because its strong suppression of 7 Be neutrinos 
helps account for the difference between the 37 CI and SuperK rates (note the x 2 values for 
the 37 C1 data in Table 0). 



5 Solutions with non-standard contributions from hep neutrinos 

Recently it has been speculated that the rise in the SuperK E e spectrum at higher energies 
could be due to a larger than expected hep neutrino flux contribution ]T5L While the 
maximum energy of the 8 B neutrinos is about 15 MeV, the hep neutrinos have maximum 
energy of 18.8 MeV. In the SSM the total flux of 8 B neutrinos is about 2000 times that of the 
hep neutrinos || , and the hep contribution to the SuperK experiment is negligible. However, 
there is a large uncertainty in the low energy cross section for the reaction 3 He + p — > 4 
He + e + + v e in which the hep neutrinos are produced, and an hep flux much larger than the 



SSM value may not be unreasonable |L5|, |lq| . A large enhancement of the hep contribution 
could in principle account for the rise in the E e spectrum at higher energies seen by SuperK. 
The hep flux normalization 7 can be determined once there are sufficient events in the region 
E e > 15 MeV. SuperK measurements of the E e spectrum in the range 17-25 MeV already 
place the upper limit 7 < 8 at 90% C.L. 0, assuming no oscillations. 

It should first be noted that although an enhanced hep contribution without neutrino 
oscillations can provide a good fit to the SuperK data for 5.5 < E e < 14 MeV, it cannot 
also account for the 37 C1 and 71 Ga rates even with arbitrary 8 B, hep, and 7 Be neutrino flux 
normalizations. The overall best fit in this case occurs with no 7 Be contribution, and has 
X 2 1 DOE = 53.7/26, which is excluded at 99.9% C.L. The contributions of the pp neutrinos, 
plus the reduced amount of 8 B neutrinos needed to explain SuperK, give a rate for the 
radiochemical experiments that is too high, even when the 7 Be contribution is ignored. 

One can also ask what happens if an enhanced hep contribution is combined with neutrino 
oscillations [IIH], although the motivation for enhancing the hep neutrino flux is not strong 
here since vacuum oscillations already can explain the rise in the SuperK E e spectrum. 
Table ^ shows the best fits in each region when an arbitrary 7 is allowed. The overall best fit 
parameters in this case are again in region C with x 2 /DOE = 32.8/25, which corresponds to 
a goodness-of-fit of 14%. The fits for most of the regions are not significantly improved from 
the standard hep flux case, and the fitted oscillation parameters are little changed. The only 
exceptions are solutions E and Z, which originally could not explain the SuperK spectrum 
as well, but with the addition of extra hep neutrinos can now also provide a reasonable fit 
to all of the data. Regions A, C, D, E, and Z are all allowed at 95% C.L. when an hep 
enhancement is included. However, the preferred values of 7 exceed the current bound from 
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SuperK, so that the role of an hep flux enhancement appears to be minimal. 



6 Summary and discussion 

The latest SuperK solar neutrino data suggest there is a seasonal variation in the solar 
neutrino flux. The hypothesis that solar v e undergo vacuum oscillations to an active neutrino 
species provides a consistent explanation of all the solar data, with a best fit given by 
oscillation parameters 5m? = 4.42 x 10 -10 eV 2 and sin 2 26 = 0.93. Oscillations to sterile 
neutrinos are ruled out at 97.7% C.L., and fits with an enhanced hep neutrino flux do not 
significantly alter the fit results. 

The existence of vacuum neutrino oscillations can be confirmed with more data from 
SuperK on the seasonal variation of the 8 B neutrino flux. The spectrum and seasonal varia- 
tions of 8 B neutrinos can also be measured in the Sudbury Neutrino Observatory (SNO) ]I7| 
and ICARUS [18] experiments. The line spectrum of the 7 Be neutrinos gives larger seasonal 



variations than 8 B ||7|, || and these may be observable with increased statistics in 71 Ga exper- 
iments, or in the BOREXINO experiment [O], for which 7 Be neutrinos provide the dominant 



signal. Accurate measurements of the seasonal variation in these experiments should be 
able to distinguish between the different vacuum oscillation scenarios || , providing a unique 
solution to the solar neutrino problem. 
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Table 1: Best fit parameters and x 2 f° r various vacuum oscillation solutions with v e oscilla- 
tions to an active neutrino. 





Parameters 








x 2 




x 2 


Solution 


5m 2 


sin 2 29 


(3 


37 C1 


71 Ga 


spectrum 


seasonal 


total 




(10- 10 eV 2 ) 
















A 


0.65 


0.70 


0.95 


4.7 


1.2 


25.8 


6.7 


38.4 


B 


2.49 


0.80 


0.79 


8.9 


1.1 


30.4 


5.4 


45.8 


C 


4.42 


0.93 


0.78 


6.7 


4.4 


17.4 


5.3 


33.8 


D 


6.44 


1.00 


0.80 


7.5 


5.2 


17.6 


6.4 


36.7 


E 


8.61 


1.00 


0.81 


7.7 


5.3 


20.4 


8.7 


42.1 


Z 


0.06 


1.00 


0.47 


5.0 


2.5 


24.8 


13.3 


45.6 



Table 2: Best fit parameters and x 2 f° r various vacuum oscillation solutions with u e oscilla- 
tions to a sterile neutrino. 





Parameters 








x 2 




x 2 


Solution 


Sm 2 


sin 2 20 





37 C1 


71 Ga 


spectrum 


seasonal 


total 




(10- 10 eV 2 ) 
















A 


0.64 


0.62 


0.96 


18.9 


1.4 


29.3 


6.7 


56.3 


B 


2.49 


0.73 


0.82 


21.2 


1.4 


33.1 


5.3 


61.0 


C 


4.42 


0.89 


0.84 


15.8 


2.7 


18.5 


5.4 


42.4 


D 


6.42 


0.97 


0.88 


15.3 


5.2 


17.3 


5.9 


43.7 


E 


8.62 


1.00 


0.91 


15.7 


5.5 


19.2 


8.1 


48.5 


Z 


0.06 


1.00 


0.48 


6.2 


2.4 


25.5 


11.5 


45.6 



Table 3: Best fit parameters and x 2 f° r various vacuum oscillation solutions with v e oscilla- 
tions to an active neutrino and a non-standard hep contribution 7 times the SSM prediction. 







Parameters 










x 2 




x 2 


Solution 


Sm 2 


sin 2 29 


P 


7 


37 C1 


71 Ga 


spectrum 


seasonal 


total 




(10- 10 eV 2 ) 
















A 


0.65 


0.70 


0.94 


7 


4.4 


1.3 


26.0 


6.7 


38.3 


B 


2.50 


0.76 


0.74 


38 


8.5 


1.1 


27.9 


5.8 


43.2 


C 


4.42 


0.90 


0.74 


35 


6.2 


3.4 


17.6 


5.6 


32.8 


D 


6.44 


0.98 


0.76 


50 


6.2 


4.5 


17.4 


7.1 


33.2 


E 


8.61 


0.97 


0.75 


66 


5.9 


4.3 


18.4 


10.3 


38.9 


Z 


0.06 


0.98 


0.44 


52 


2.7 


2.2 


19.3 


14.9 


39.1 



7 



Cl+Ga Cl+Ga+SK spectrum Cl+Ga+SK spectrum 

SK spectrum SK seasonal +SK seasonal 




0.01 I 1 1 1 1 ] I 1 1 1 1 I I 1 1 1 1 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

sin 2 (20) sin 2 (20) sin 2 (29) 



Figure 1: Regions in sin 2 29-5m 2 parameter space (a) allowed at 95% C.L. by the 37 C1 
and 71 Ga capture rates (solid curve) and the Super-Kamiokande E e spectrum (dotted), (b) 
allowed at 95% C.L. by the combined 37 C1 and 71 Ga capture rates and Super-Kamiokande 
E e spectrum (solid) and excluded at 68% C.L. by the Super-Kamiokande seasonal data 
(shaded), and (c) allowed at 95% C.L. by all the data. 
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Figure 2: Ratio of predicted E e spectrum in Super-Kamiokande to the SSM for (a) models 
favored by a global fit, C (solid curve), D (dashed), and A (dotted), and (b) models disfavored 
by a global fit, B (solid), E (dashed) and Z (dotted). 
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Figure 3: Ratio of predicted rate to SSM prediction versus time of year in Super- 
Kamiokande for (a) models favored by a global fit, C (solid curve), D (dashed), and A 
(dotted), and (b) models disfavored a by global fit, B (solid), E(dashed), and Z (dashed). 
Results for no oscillations with 8 B normalization of 0.47 are shown by the dash-dotted curves. 
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